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The two dimensional (1+1) sine-Gordon model finds many 
applications in condensed matter physics. These in turn pro- 
vide an experimental means for the study of topological de- 
fects, some of which may have had a huge impact on the early 
universe. As a first step in trying to exploit this analogy and 
also others which exist with low-energy QCD, we study boson- 
isation in the massive Thirring and sine-Gordon models at 
finite temperature T and nonzero fermion chemical potential 
fj,. Both canonical operator and path integral approaches are 
used to prove the equality of the partition functions of the two 
models at T > and = 0, as was recently shown. This en- 
ables the relationship between thermal normal ordering and 
path-integral renormalisation to be specified. Furthermore, 
we prove that thermal averages of zero-charge operators can 
be identified as long as one uses the usual T = identifica- 
tion between coupling constants. Analysis of the point-split 
regularised fermion current then leads to the thermal equiv- 
alence between sine-Gordon kinks and Thirring fermions. At 
^ 7^ and T > 0, we show, in perturbation theory around 
the massless Thirring model, that the bosonised theory is the 
sine-Gordon model plus an additional topological term which 
accounts for the existence of net fermion charge excitations 
(the fermions or the kinks) in the thermal bath. This result 
generalises that recently obtained for the massless case, and it 
is the two-dimensional version of the low-energy QCD chiral 
Lagrangian at finite baryon density. 



I. INTRODUCTION 

The sine-Gordon (SG) theory in two dimensional 
Minkowski space-time with metric (+, — ) is described by 
the Lagrangian density 
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C-sgW\ = ^df,(j)d^'(t) + ^ cos A(/) - 70, 



(1) 



*Proceedings of the talk based on [1] and given by D. A. Steer 
at the 5th International Workshop on Thermal Field Theory, 
Regensburg, Germany, August 10-14, 1998. 

^ Gomez@eucmax.sim.ucm.es 

*R. Rivers@ic.ac.uk 

^ D.A.Steer@damtp.cam.ac.uk 



where ^ is a real scalar field, 70, ao and A are bare 
parameters to be renormaHsed later, and = {t,x), 
H = 0,1. Notice that apart from the usual kinetic term, 
the potential term is periodic so that there are an infinite 
number of degenerate vacua whose value, depends on 
the coupling constant A: (j)^ = 2n7r/A where n a 1L. 
We comment that the Lagrangian is invariant under 
— > (j) + (j)y and <j) ^ —(j) meaning that the sign of ao 
is unimportant, and the name sine-Gordon is related to 
the fact that the equations of motion contain a sin <f> term 
(see il). 

The massive Thirring (MT) model, on the other hand, 
is a model with a fermionic field ip = {ipi, V'2) and a four 
Fermi interaction with coupling g: 

CTh[i^,iP]^i^{P-mo)^^ + ]^gj^,{x)f{x), (2) 

where jfj,{x) — ip{x)^fj^ip{x) and mo is a bare mass. For 
positive coupling constant, g > 0, the interaction term 
is attractive and there are fermion- antifermion bound 
states. 

Though it maybe somewhat surprising since one is a 
model with a bosonic field and the other one with a 
fermionic field, it is very well known that these two the- 
ories are linked. In particular, they provide one of the 
earliest examples of duality in which the weak limit of 
one theory describes the same physics as the strong limit 
of the other and conversely since 



A^ 

47r 
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l-fg/Tr' 



(3) 



Hence perturbative calculations in one theory tell us 
about non-perturbative effects in the other. Identity (H) 
was proved in perturbation theory about ao = = too 
(see below) at T = u = by Coleman using canonical 
operator methods ||^ and later on using path integrals 
1^,^. Note the special value g — Q in which the Thirring 
model is a free theory of massive fermions, corresponding 
to A^ = 47r in the SG model. Also the MT bound states 
exist for 5 > ^ A^ < An. 

The duality (^) is directly related to the fact that 
the models exhibit bosonisation, in which a theory of 
fermions is equivalent to a theory of bosons. In general 
the bosonic (fermionic) theory may also have fermionic 
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(bosonic) excitations; this indeed occurs in the SG and 
MT models and we will see another example below in 
low-energy QCD. In the case of the SG and MT models, 
bosonisation schematically works in the following way 
(see for a summary as well as a list of the relevant 
articles). Consider first the SG model. Perturbations in 
A about one of the minima of the degenerate potential 
give rise to the usual bosonic simple harmonic oscillator 
spectrum. The model also has other excitations. Recall 
that there are topologically non-trivial solitonic solutions 
to the classical equations of motion, and that these arise 
from the vacuum degeneracy: there is no reason why at 
x — > cxD the system should be in the same vacuum as at 
X —oo. As a result one can have finite energy solutions 
with non-zero charge 

/oo \ c) A 

(4) 

The simplest 'kink' solution has Q = 1 and an energy 
proportional to ao/A^, whilst the anti-kink (Q = —1) is 
obtained from the kink by taking (j) — > —0. Furthermore, 
since the system is integrable, exact expressions for multi- 
kink and anti-kink states are known and one can con- 
struct moving kinks by Lorentz transformation. These 
may then be quantised using semi-classical methods . 
There are also classical solutions with Q = corre- 
sponding to 'breather' solutions in which kinks and anti- 
kinks oscillate periodically about each other j^. These 
can be quantised by WKB methods The relation- 

ship between these SG excitations and the fundamental 
fermionic excitations and bound state excitations of the 
MT model is summarised schematically in table I. 



MT MODEL 


SG MODEL 


fermionic fields 


bosonic fields 


9 


A ^ 1/g 


ff bound states 


bosons, kink - anti-kink breather 


(5>0) 


solutions (A^ < in) 


/ 


kinks 




™o 7^ ^ see later 




mo = ^ topological term |^ 



TABLE L Table showing schematically some of the 
links between the SG and MT models at T = 



We see that the kinks themselves are fermion-like excita- 
tions corresponding to the fundamental fcrmions of the 
Thirring model. For A^ < 47r, the bosonic {Q = 0) exci- 
tations of the SG model correspond to the bound states 
of the MT model. This is what is meant by bosonisation. 

What is the motivation for studying these models at 
T > and /Lt ^ 0? One first reason is that these mod- 
els could be used to develop techniques which may then 
be applicable to more realistic 4D theories, in particu- 
lar QCD. Recall that at low temperatures, where the 



quarks and gluons are strongly confined into hadrons, one 
can describe the system by an effective chiral bosonic la- 
grangian (CBL) for the lightest mesons (pions, kaons and 
eta) which are the Nambu Goldstone bosons (NGB) of 
the chiral symmetry breaking. The relationship between 
this chiral bosonic theory and the original fermionic QCD 
has many similarities with the relationship between the 
SG and MT models (see table II). 



QCD 


CBL 


quarks and gluons 


bosonic fields 


Us (large at low 


small expansion parameter 


energies) 


Ax ~ 47r/^ ~ 1 GeV 


qq bound states 


TT, K, 7] (NGB of chiral 
symmetry) 


baryon 


skyrmion (topological defect) 


+ fi for baryons 


topological term j|] 



TABLE II. Table showing schematically some of the 
links between QCD and the low energy chiral 
lagrangian. 



For example, the bound states of QCD are the mesons 
of the chiral theory whilst the baryons of QCD corre- 
spond to skyrmions — topological defects [llO|jll[ — in the 
chiral theory which to lowest order is the non-linear sigma 
model ijl^ (c./. the relationship between the SG and MT 
models in table I). Note also that the SG Lagrangian in 
(|l|) corresponds to a non-linear sigma model in 2D for a 
single Nambu-Goldstone-like field </>. Although there is 
no spontaneous symmetry breaking in 2D [p^ , the poten- 
tial term in (|l|) breaks explicitly the symmetry <j) ^ <j) + a 
with a G IR (which we will call the chiral symmetry for 
reasons to become clear below) still preserving the sym- 
metry (j) (f) + (j)y. These two symmetries are, respec- 
tively, the counterparts of the chiral and isospin sym- 
metries for QCD, ao and A playing the role of the pion 
mass squared and the inverse of the pion decay constant 
/tt respectively. On the other hand, the chiral symme- 
try transformations in terms of the Thirring fermion are 
-0 — s- exp(j 075 )■)/'• The massless Thirring model is chiral 
invariant, the fermion mass term breaking that symmetry 
in the same way as the a term does in the SG Lagrangian. 

A second motivation for studying these models at 
T > and /i 7^ comes from cosmology. The SG 
kinks are the 2D analogues of cosmic strings, line-like 
defects formed in 4D when a system goes through a sym- 
metry breaking phase transition (say of some group G 
to a subgroup H) for which the first homotopy group 
TTi{G/H) ^ 0,111 . The strings trap regions of the un- 
broken high energy phase and so have energies per unit 
length which can be very large, depending on the critical 
temperature. In the context of the early universe and 
cosmology, cosmic strings may have played an important 
role because as the universe cooled and expanded after 
the big bang, it went through a number of phase tran- 
sitions some of which may have led to the formation of 
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strings. In particular, strings formed at the GUT phase 
transition have huge energies per unit length {n ~ 10'^^ 
GeV^) and hence significant gravitational efi'ects, and so 
it is thought that they may have been responsible for the 
temperature fluctuations in the cosmic microwave back- 
ground radiation and for seeding the perturbations which 
led to the formation of structures such as galaxies Jist . 

However, in order to make detailed predictions as to 
their effects, it is important to know the initial distri- 
bution of strings and whether or not it contains infinite 
strings. This is a very difficult task to undertake analyt- 
ically , but recently progress has been made by using 
the analogy between experimentally observable systems 
such as '^He and '''He and the early universe [0. It would 
seem, however, that there is an even more simple ex- 
perimentally accessible system with which one could try 
to test ideas of defect formation, and that is a Joseph- 
son junction This device consists of two layers of 
superconductors separated by a thin dielectric barrier, 
typically of the order of 5nm. Denoting the macro- 
scopic wave function of each of the superconductors by 
= l^al exp{i6'Q} where a = 1,2 labels each of the 
two layers, then Josephson tunnelling of the Cooper pairs 
across the dielectric layer results in A9 = 9i — 62 satisfy- 
ing a SG equation ||l8|. As the two superconductors are 
taken through the phase transition, kinks are formed in 
the junction and these are observed experimentally p9| . 
Experiments are also done to see what is, for example, 
the effect of the geometry of the set up on the kinks and 
their dynamics ||20|| and these devices are used as sources 
of some of the highest energy micro-waves finding ap- 
plications in satellites |lj,|l^. One idea might therefore 
be to see whether one can indeed use such experiments 
with Josephson junctions to test ideas of defect forma- 
tion. The hope is in particular that the situation can be 
much simplified in this case because of the duality be- 
tween the SG and MT models. First, however, one has 
to check what is the relation between these two models 
at r > and fJ. ^ 0- In particular, could it be possible 
that the relationship between coupling constants (^) is 
a function of temperature? As we will comment in the 
conclusions, these models also provide examples of other 
phenomena which may too be more easy to study in 2D 
rather than 4D. 

Here we report on our first steps in these directions. In 
section || we briefly summarise some of the particulari- 
ties and basic results of 2D which will be useful to bear 
in mind for the rest of the work. Result s arc then given 
in section HI. First, in subsection III A, we outline the 
main steps which must be taken to extend Coleman's 
work Q on the SG model to T > using an operator 
approach. The results of showing the equivalence 
of the partition functions are reproduced though the ap- 
proach is entirely different. Path integral methods are 
used in the rest of the work. We show that not only the 
partition functions are equivalent but also that thermal 
averages of correlators of zero-charge operators evaluated 
at different space-time points coincide. The relationship 



between normal ordering and regularisation is specifled 
and we analyse the T > situation in which there is a 
net number of fermions, ^ ^ 0. The partition function 
is calculated in that case so extending [|j to the mas- 
sive case and to ^ 7^ 0. The analogy between these 
models and a classical gas of particles is then noted in 
the conclusions. Details of the calculations and results 
presented here may be found in references [0J2|]. 



II. PECULIARITIES OF 2D 

The following three basic, though important, results 
hold for free fields in 2D at T > 0. They will be useful 
in the following sections. 

1) Equality of free massless boson and fermion 
partition functions 

Consider a free bosonic field of mass /2 and a massless 
free fermionic field. The partition functions can be ob- 
tained by writing the Hamiltonians in terms of a normal 
ordered part plus an infinite vacuum energy; they are 
respectively 



\nZi;^^{T) = -L I ^ 



-j-oo 



dk 



■ + 



J-oo 27r 



ln(l - e-^"^"-^^] 
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In 1 



(5) 



(6) 



Here ujk 



k'^ + p,'^, /? is the inverse temperature f3 — 



1/T, and L is the spatial dimension which we are taking 
to infinity. As usual, all the thermodynamic observables 
are obtained from the logarithm of the partition function 
and its derivatives by dividing by f3L. 
Observe that when jl = 0, 



^Inf 
27r V 



1 - e 



^-2rflnfl + e-^NU-^, 



2tt 
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and hence it follows that ignoring vacuum terms, the two 
partition functions (||)-(||) are equivalent for jl = 0: 



exp 



TiLT 
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(7) 



Thus one cannot tell the difference between the bulk 
quantities of these two systems: equality is perhaps 
the simplest example of bosonisation. As will be seen be- 
low, the equality of the SG and MT partition functions 
at T > rests on (|^) since in those models we work in 
perturbation theory, expanding about ao = in the SG 
model and about mo = in the MT model (and hence 
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about massless bosonic and fermionic theories). 
2) Free thermal propagators 

We will need the free propagators for boson and fcrmion 
fields at T > 0. Again working with a boson field of 
mass jl and a massless fermion field, and calculating in 
the imaginary time formalism with i = it, these are found 
to be respectively pif 



At{x) 

Sai3{x) 



1 

47r 
f 



2(3 Q^x) 



(8) 
(9) 



where we have expanded the boson propagator about ji 
0, and X is a constant. The Q variable is given by 
Qo + Qi where 



Qo{x,t) = -cosh( — )sin( — ), 

— 7TX irt 

Qi{x,t) = -sinh( — )cos( — ), 



(10) 



so that Q{x) is a Lorentz scalar. The indices a, (3 are 
Dirac indices and we have worked with the 2D Euclidean 
7 matrices 



7° = 



75 



so that 



{7^7^} = 25'^'', 7A.75 = -^£^-7. 



with eoi — 1. Observe that (||) is both ultra-violet diver- 
gent {x —>■ 0) as well as infra-red divergent (x — > cxd or 
// — > 0) since 



Q\x,t) 
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(11) 



An interesting property of the propagators (||)-@ is that 
they can be directly obtained from the corresponding 
T = propagators through the substitution x^ — > /3(5^. 
This is in fact true for all contours in the complex time 
plane, and can be viewed either as a result of the con- 
formal invariance of the free theories, or as arising from 
solving the Green function equation for the Coulomb po- 
tential on a cylinder of radius (3 p^ . 

3) Ultra-violet divergences and normal ordering 

Finally, recall that the UV divergence structure of 2D 
bosonic theories is much simpler than that of 4D ones. 
The reason follows from the fact that in d dimensions 



with an interaction of the form (j)^ , a diagram with n ver- 
tices and E external lines has a UV degree of divergence 
L> of 



D 

So with d = 2, 



d 

--l]E + n 



-{d-2)-d 



D = 2-2n 



and the only divergent diagrams Vr are tadpole dia- 
grams. Since the SG lagrangian contains a term cos = 
^(— l)^(/)^"/(2n!) then in the operator formalism it fol- 
lows that all UV divergences of the theory should be re- 
moved through normal ordering (as this removes tadpole 
diagrams). In path integral methods there is no operator 
normal ordering, and the divergences are removed using 
different methods (see below). 



III. RESULTS 

A. Bosonisation in the canonical operator approach 

at r > and = 

We have used operator methods to extend the paper 
of Coleman to T > and fi = 0. This enabled us to 
prove that the partition functions of the SG and MT 
models are identical 



ZsG{T,n^O) = ZMTiT,n = 0) 



(12) 



provided the coupling constants of the theories satisfy 
(^) which is temperature independent. In fact, ( [l^ ) was 
already proved using ve ry diff erent path integral methods 
in 1 21 1 (see also section [II B ). 

In the SG model the main steps in the calculation 
Zsg{T, M = 0) ^-re the following: 

1 . Remove all U V divergences by normal ordering the 
SG Hamiltonian. 

2. Deal with the IR divergence of the propagator (H). 

3. Calculate the partition function Zsg, a sum of ther- 
mal expectation values of free (interaction picture) 
operators each of which have been normal ordered 
as a result of step 1. 

We now outline the main features of each of these steps, 
and also comment on the differences between this T > 
case and the T = one discussed in jj] . 

Step 1: Removal of all UV divergences by nor- 
mal ordering 

As stated in point 3 above, the partition function con- 
tains thermal expectation values (TEV's) of normal or- 
dered operators. Therefore to simplify its calculation, we 
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remove UV divergences by using thermal normal order- 
ing (TNO) introduced in |24| rather than standard zero 
temperature normal ordering (which places annihilation 
operators to the left of creation operators). By construc- 
tion TNO, which is denoted by N^^ , guarantees that for 
any^ operator (other than the identity) • in the interac- 
tion picture, 



<iV-^'^H >o =0, 



(13) 



whereas for usual normal ordered products 
■C: • :>o 7^ 0. 

The operation N^"^[»] was defined in to place the 
"positive" parts of the operator, •+ (a combination both 
of annihilation and creation operators) , to the right of the 
"negative" part, See |2j] for the exact definitions. In 
the above equations the angular brackets denote a ther- 
mal expectation value, and the subscript zero indicates 
that we work with free (interaction picture) fields, and 
hence in perturbation theory. Identity (^3|) will greatly 
simplify calculations below. 
The SG Hamiltonian is 



HsG = 



dx 



Tr\ 1 d(j)\ ao .7 
Y+2\d^ -^cosA0-7o 



Ho 



dx cosX$ + 7o 



(14) 



which must be divided into a free and interacting part so 
as to apply TNO. Although the term cosAc^ itself con- 
tains a mass term on expansion in powers of A, we want 
to keep A of arbitrary size. Consider therefore 



HSG = 



dx [ ^ cos X(j> + + 70 



so that perturbations are about a scalar field of mass p. 
To take account of this fact, TNO is now denoted by 
Np^[»] and similarly we add an extra mass label to the 
propagators (||); At{x) AT{x;fl). Hence in the case 
of HsG above we will be dealing with At{x] p). The link 
between p and the regularisation scale in path-integral 
methods is discussed in section 



IIIB 



TNO of ( |r^ ) may be carried out by using the identity 



J^d^£c J^tfyj{x)ATix-y;p)j(y) 

where Tc means contour ordering Had we used T — 
normal ordering, the zero-temperature propagator would 
have appeared in (|l^) rather than the finite temperature 
one: TNO means that the Q variables of (|lO|) are built in 
from the start of the calculation. Finally the UV diver- 
gence of At{x; p) is regulated by cutting off the theory 
replacing 

ATix;p) At{x-p-K):^At{x-p)-At{x-K) 

(16) 

where A is a large mass. Hence the constant K in 
the propagator (||) cancels, and At{x\p\K) is now both 
non-singular as well as (3 independent for x — > 0; 
At(0; p; A) - - In (p^A^) . Combining (§]) with (|l|) 
for j{x) = \5{x — y) gives, for example 



A2 



N 



ES 



iX(j>{y) 



(17) 



After similar manipulations, (|l^) can be written as Q 



HsG - N. 



ES 



dx cos X(f) — L"f 



(18) 



where ao has been multiplicatively renormalised, 



and 7o has been renormalised according to 
7 = 7o - Et{p) 

where 

1 



(19) 



(20) 



ETip)^^{[7T+.n-] 

/■+°° dk 1 



(900)+, (900)- 
2k^ 



Here Eq{p) is an infinite temperature independent con- 
tribution; Ea{p) = L Hi ^77^- The temperature 
dependent part proportional to the Bose-Einstein distri- 
bution Nk^p = (e^'^'^'P — 1)~^ is finite. The coupling A is 
unchanged. 



j,ES 



^This could be the field operator 0, the momentum operator 
■if or any other composite operator. 



These equalities may at first sight seem surprising, but they 
might be clarified by noting that here the advanced and re- 
tarded thermal propagators are equal so that A^{x — y,p) = 
Arix — y, p) = [(f>'^ (x) , (b~{y)] where the positive and negative 
parts refer to those of 124]. 
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Thermal normal ordering the SG Hamiltonian has 
therefore absorbed all UV infinities just as zero tempera- 
ture normal ordering does Q , but it has also introduced 
some extra T-dependent finite terms. 

Steps 2 and 3: IR divergences and calculation 
of the partition functions 

The IR divergence of the boson propagator (||) is removed 



by introducing a mass fl into the SG Hamiltonian (18). 
At the end of the calculation we let jl —> 0, and hence 
are free to add the extra mass term within the normal 
ordering giving the Hamiltonian 



H = 



ES 



=:io- 



ES 



(^H(^ - jLj ~ ^ y ^cos X(j) 
/ dx(cosA0j 



Ai. 



(21) 



Here 



Hq + dx^fl^(i>^, and Aq and Aj denote 
the free (a = 0) and interacting Hamiltonians respec- 
tively. In perturbation theory the SG partition function 
is therefore given by 



ZsciT^fJ, = 0) = lim Tr 



lim > — 7 / dti . 

„ n! Jo 

{e-^^«r, [i,(ti)...i,(t„)]}, (22) 



dt„ X 



Tr i e 



so that one only needs calculate free expectation values 
< • >o= Trexp{-/3Ao»l_/Trexp{-/3lo}-^ Indeed, the 
correlator appearing in ( |22|) can be obtained by showing 
first that lEIl 



ES 
P 



>0 



ES 



>0 X 



xY[[(3^fi^\Qix,-Xk)\' 

j>k 



In fact one has to be extremely careful when considering the 
precise form of the thermal weight exp{~/3Ao}. The reason 
is that it contains two different mass scales p and jl, meaning 
that Ao is not obviously diagonal and so not obviously of the 
form J dkio(k)a^ {k)a(k) as is usually the case and as was 
assumed in 124]. See M for details about this point. 



sir X]j 



nn[/3viQ(^ 

j=i j>k 



Xk) 



(23) 



where we have expanded the exponential inside the nor- 
mal ordered term and noted that all terms which contain 
vanish (by (|l3|)) apart from that with n — 0. Then 
observe that the terms proportional to jl in ( p3| ) have 

a contribution Thus in the limit /i — > 0, 

only configurations with "^Xj =0 contribute (a condi- 
tion which will become analogous to the fermion chiral 
selection rule of section pIB ). In (^3|) n must therefore 
be even. If we let m 
and define 



n/2, yj — Xj for j = n/2 + 1, . . . , n 



A+ = N, 



ES 



(24) 



then from (H) 

■m 

^TcY[A+{xj)A^{yj) >o 



n 



21 4^ 



nr=i[/3viQ(^,-2/fc)p]^ 

Hence the partition function (|2^) is |Q] 

2 



(25) 



1 



n 

n">fc [IQi^j " Xk)\^\Q{yj - yk)\ 



a 

2A2 



n 



21 t;; 



d^Xjd^yj X 



(26) 



where Jrpd^x = Jq dt 

Trje-^'M = e'^T«^Z°'^(r) which is finite for 70 sat- 
isfying (!§ and < 47rJ^ 



For the MT model, calculation of Zmt(T,ij, = 0) in 
the operator approach is rather more complicated (it is 
based on a paper by Klaiber [^) and is discussed in [0. 
In perturbation theory one can show that the partition 
function ZMT{T,fi = 0) is given by [Q 



ZMT{T,tl = 0) = 



n=0 



2/3 \t) 



2n 



Ilk<j [Q^i^j - xk)Q^{y3 - Vk)] 



riLi [Q'^ixj-vk. 



1 l — K^ I TT 



n / d^^jd'vo X 
i=i ■'^ 

(27) 



*From the behaviour of the Q variables in (|ll|), one can see 
that for y? > 47r there are extra divergences in (^6|). The 
treatment of these is commented on in |lj . 
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where we have chosen to renormalise the MT model at 
the scale p. Here is given in 

(28) 



and m is the renormalised mass (see and also section 

inBl) 



m = mo (A/p) (29) 

Thus term by term ( p6| ) and (|7|) are identical provided 
that a) the parameters of the two theories are identified 
as in (0), b) that a and m are related by 



(30) 



and c) that 70 = as then [T) = Z^ (T) by (|). 

We have therefore extended the work of Coleman 
to T > 0, /i = using operator methods. This was ren- 
dered more simple through the use of TNO. The results 
obtained are in fact the same as those of |^ though the 
approach has been entirely different. We now turn to 
path integral methods. 



B. Path Integral bosonisation at T > and p. = Q 

We have used path integral methods not only to prove 
the equality of the two partition functions as in ( [1^ ) 
but also to prove the equivalence of certain sets 
of correlators of operators evaluated at different space 
time points. These correlators cannot be obtained from 
the partition function which only contains information 
about global thermodynamic observables like the pres- 
sure or the condensates, but not about correlators which 
physically yield, for instance, thermal correlation lengths. 
We outline the main points of such PI calculations for the 
SG model and then the MT model. This will enable the 
link between normal ordering and PI regularisation to be 
made. Our results for the correlators will also be stated 
more precisely. 



1. SG model 

As always in path integral methods, one works with 
the generating functional. Once again it is useful to start 
with free boson fields and in particular to calculate the 
correlator (|2^) which may be obtained from the free bo- 
son Euclidean generating functional 

1 



Z^[J-T]^N0 I a<p exp. 

'periodic 



cPx^ X 



J{x)AT{x-y)J{y)}. (31) 



Here is an infinite T-dependent constant arising in the 
path integral description |26|, the propagator is given in 
(^ for small p(3 and the free boson partition function is 
Z^{T) = Z^[0;T] as in @. Note that we have removed 
the p labels on propagator and partition function as there 
is no longer any possible confusion with other mass scales. 
Now define 

(remember that we do not normal order the operators in 
path integral and so this differs from the definition (p4|)). 
Then 



< [| A+{xj)A.{yj) >o = 



ZiiT) 



with 



n 

J{z) = -^\J2[s^^\z-x,)-S(^'Hz-y,) 



As opposed to section |III A| , here the UV divergence of 
the propagator (@) is regulated by replacing 



g^(0, 0) ^ g^(eo, ei) = -t- O(e^) 



(32) 



where 0+ and = n'^ {e^ + ef) . From (^, (||) and 

(Mj it follows that 



n,>fe [Q'^i^j - Xk)Q'\yj ~ yk)] 



A^/47r 



nLi [Q'{x,~yk)] 



A2/4ir 



(33) 



The correlator is divergent due to the short-distance 
(UV) divergent behaviour of the composite operator 
exp[z(/)(a;)], which needs to be renormalised. We do this 
in the usual way through the replacement of exp[ia(/)(a;)], 
with a G IR arbitrary, by 



[exp [ia0(x)]] 



hare 



{ep) 



exp [ia0(x)]] 



(34) 



Here p is an arbitrary renormalisation scale and the su- 
perscript R will denote renormalised operators. Note 
that (34) is analogous to (pT|) with the identification 



1 

^=A 



(35) 



though in the operator formalism the renormalisation 
was carried out through TNO. Also observe that the 
p's appear in the same way though they have different 
origins — in the operator approach p corresponded to an 
arbitrary mass at which normal ordering was performed 
whereas here it is the arbitrary renormalisation scale. 
From equations (|3^) and (|3^) observe that (|3^) reduces 
to (p5|) as required. 
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In the full SG model one again works with the gener- 
ating functional which is expanded formally in powers of 
oio/}? The partition function Zsg{T,h = 0) is just 
obtained by setting the external sources to zero. We find 
[Ql that with the regularisation ( |3^ ) of the propagator, 
ao is renormalised just as in (|9|) for all the divergences 
to be eliminated. Once again the partition function is 
given by (H). 



the sources are set to zero. Its calculation requires some 
standard manipulations (|,^,^,^ (writing the quartic 
Thirring interaction as a 'gauge-like' interaction as well 
as calculating the axial anomaly |^^) and one again ob- 
tains the result (p7|). 

3. Zero-charge operators equivalences 



2. MT model 

In calculating the MT partition function in perturba- 
tion theory about mg — 0, the correlator analogous to 
( p5| ) in the SG model is just the TEV of insertions of the 
operators (J±{x) = ip{x)P±'4>{x). Note, however, that the 
massless fermion theory is invariant under chiral trans- 
formations '0 exp(ia75)'0- Under such transformation 
u±{x) exp(±2iQ!)(T-|-(2;) and therefore the thermal av- 
erage of a product of a± (x) operators will vanish in the 
massless case unless the number of cr+ and ct_ is the 
same. This is the chiral selection rule, which only holds 
for mo = 0. Following ||2^], the required correlator is 
obtained by shifting tj) —>■ ipj'^ so that, naming ipa with 
a = 1, 2 the two components of the bispinor -0, the free 
massless theory decouples into two free theories for the 
spinors ijja, and we have cr+ 4'2'4'i and a- — > V'i'02- 
One obtains [0j27) 



i=i 



-Vk)] 



(36) 



Notice that the above correlator has exactly the same 
structure as the boson correlator (p3) — this is another 
peculiarity of 2D. However, unlike (|33|), (|36| ) is finite since 
it contains no product of fields at the same space-time 
point and there are no mixing terms between ipi and ip2 
in the Lagrangian. 

For the MT model one does have to worry about renor- 
malisation (section III A). The reason is that whilst the 
chiral symmetry is still unbroken so that (t±{x) correla- 
tors still appear in the same number if g 7^ and mg = 0, 
now (ipjf^ip)'^ — > 4:ipiipiip2^2 when — > ipj'^ and there- 
fore there is mixing between ipi and "02. Thus products of 
fields at the same point appear and the (j± (x) correlator 
becomes divergent: in the same way as the boson op- 
erator exp[ia0(x)], the a±{x) composite operators need 
renormalisation. Also as in the SG model, those are the 
only infinities we have to worry about and they are ab- 
sorbed in the renormalised mass m as in (^^ whilst 



(37) 



Given these renormalisations, the MT model partition 
function is obtained from the generating functional when 



So far we have shown the equivalence of the SG and 
MT partition functions with n = using both operator 
and PI techniques. Further equalities between correlators 
of different operators in each of the two models may also 
be shown to hold at T > and fi = — we simply state 
the results here. Further calculational details may be 
found in Q]. 

The equivalence of (^ ) and (^6|) in the free massless 
bosonic and fermionic theories (after renormalisation) 
may be extended to the SG and MT models where it 
becomes 

N 



2N ^ 

(0 «Tcn^?(^.M''(%)»5G. 



We have also considered more complicated cases^. There 
the calculation goes through in a similar way and leads 
to the expected result in which 



(38) 



When dealing with correlators including insertions of 
the current operator (x) in the MT model one is again 
faced with a product of field operators at the same point, 
and hence with additional divergences. Using point- 
splitting regularisation ||2^ and taking particular care 
to ensure that the Ward identities are satisfied, we have 
proved that there are no extra divergences to renormalise 
and 

< T,jj^-{x)a^{zi)a^{z2) ~>mt = 

T^e^^ « T,aXa:)A«(zi)A«(^2) »5G 

where j^ (a;) is the regularised current jl]. More generally 
we find that the identities which hold between different 
correlators seem to be simply obtainable through the re- 
placements (Bq) and 



^For example there could be unequal numbers of a+ and 
CT_'s, since the chiral symmetry is broken in the MT model 
with mo 7^ 0. 
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J 



ZTT 



(39) 



which are usually called operator bosonisation identities. 
We stress, however, that firstly we do not believe these 
identities to be strong — that is, we only expect them to 
hold between thermal expectation values — and secondly 
it is not immediately obvious that they hold inside all ex- 
pectation values, though it does certainly seem to be the 
case for the zero fermion charge ones we have considered. 



C. Bosonisation of the massive Thirring model at 

T > 0, ^ / 

Finally we have studied the Thirring model at nonzero 
chemical potential /i for the conserved charge Qp — 
J dxj^ (i, x) which is the net number of fermions minus 
anti-fermions; 

CTh[i',t,fi] = -^{^ + mo)^P + ^g^jaix)rix) + fif{x), 

and we have calculated the grand-canonical ensemble 
partition function Recall that now the averaged net 
fermion number density p(^) = (/3L)~^ <^ Qf ^ is no 
longer necessarily zero, and so a natural question to ask is 
what is the bosonised version of this theory. The answer 
was obtained in |^ for the massless case where the free 
boson partition function acquires an extra /i-dependent 
term. In the massive case with partition function 



dijjdip exp 

antiper 



we find in perturbation theory around the massless case, 
after using some of the results in |^ and performing some 
functional integral manipulations, that 



ZTh{T,n) = Z5G^(T,^) 



where 



d<p exp 

periodic 



(fx {CsG 



(See for a discussion of the boundary conditions.) 
That is, the bosonised action is the SG model plus an 
extra term which is topological in that it only depends on 
the value of the field at the spatial boundary {x = ±oo). 
From (^) this term is interpreted as the result of excita- 
tions with net kink (fermion) charge being present in the 
thermal bath and having associated chemical potential fi 
in the grand-canonical ensemble. Recall that an analo- 
gous contribution was found in the chiral Lagrangian for 
low-energy QCD where the role of kinks (Thirring 
fermions) is played by the skyrmions (QCD baryons): 
the chiral Lagrangian at finite baryon density acquires, 
amongst other things, a new factor fiQsK with Qsk the 
skyrmion topological charge. 



IV. CONCLUSIONS AND OUTLOOK 

After motivating the study of the SG and MT models 
at T > and /i 7^ 0, we have summarised some of the 
results obtained in . Firstly, with zero chemical po- 
tential we were able to show using both operator methods 
and path integral methods that the partition functions 
were equivalent. This also enabled a link to be made be- 
tween the arbitrary scale p at which we carried out ther- 
mal normal ordering in the operator approach and the ar- 
bitrary renormalisation scale in the path integral — these 
are identical. We then studied correlators of operators at 
different space time points in each of the models. Such 
results will be of crucial importance for the application of 
this work to the estimation of the number of topological 
(kink) defects formed in a phase transition, an extension 
of this work which we motivated in the introduction [ p2[ . 
Also of relevance to this future work is the relationship 
between the two models in the presence of a net n umber 
of excitations; we stated our results in sections IIIB3 



and III C 



Finally, we are studying the link between the sine- 
Gordon model and a ID classical gas of positive and neg- 
ative charges with non-zero fugacity [g2[. As a result of 
this (highly studied system) and of the relationship be- 
tween the sine- Gordon model, Josephson junctions and 
the massive Thirring model, we would tentatively sug- 
gest that such models should perhaps not be overlooked 
as ones in which to develop or test calculational methods. 
For example, one could try to investigate such difficult 
quantities as the pressure or the fermion number den- 
sity, which could then give some insight into real physi- 
cal problems. Similarly we hope to investigate the precise 
nature of the transition at — Stt. 
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